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ABSTRACT 



A harmonic analysis of the observed diurnal variation of the 
coefficient of eddy viscosity is conducted at selected levels followed 
by a mathematical analysis of the amplitude coefficients in order to 
determine a function of height and time needed to represent this 
variation,, It is found by this approach that a lengthy series of terms, 
consisting of exponential and polynomial functions of height and 
trigonometric functions of time and height, would be needed to represent 
the observed variation, 

A two layer empirical model is developed from the observed data 
which gives the coefficient of eddy viscosity as a function of stability 
and height in the lower layer, and as a function of vertical shear of the 
horizontal wind and height in the upper layer. 

The author wishes to express his appreciation for the assistance and 
encouragement of Professors G„ J„ Haltiner and A, B„ Mewborn of the 
Uo S„ Naval Postgraduate School in this investigation. 
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Introduction, 



In an attempt to account for the observed wind structure in the 
friction layer, meteorologists have devoted considerable effort to 
determine the variation of the coefficient of eddy viscosity with 
respect to time and space. 

The layer of frictional influence has been subdivided into the 
laminar, surface, and spiral layers each possessing certain 
characteristic properties c The coefficient of eddy viscosity, , 

has been shown to vary linearly with height in the neutral surface 
layer and as a power function of height in the non-neutral surface 
layer. Using wind data obtained from double theodolite observations 
and assuming a steady state, Mildner |^6^j solved the horizontal 
equation of motion for J ' utilizing finite differences for the 
derivatives of velocity 0 He showed that /A increased up to a fairly 
low elevation in the spiral layer (2^0 meters) and then decreased to 
a small residual value near the gradient level. 

The Taylor [~ 7~1 solution for the wind in the spiral layer assumes 
A constant, along with a number of other broad assumptions, and 
yields a general representation of the observed vertical wind variation. 
More recently, Buajitti and Blackadar j 2*] and Haltiner |_ 4 J have 
sought the distribution of eddy viscosity as a function of time and 
height in order to account for the diurnal variations of the wind 
structure in the friction layer c Haltiner [4j| assumes the 
coefficient of eddy viscosity to have the form 

K s JA = g(t) k(z) 



1 



v/here 



g(t) - 1 + b (SinOJ t + c Sin 26J t) 

and k(z) takes on several different exponential foras 0 

It is the purpose of this investigation to determine the degree 
of complexity of the function of yM (z, t) needed to yield observed 
values of yM ; and also to determine, if possible, the coefficient 
of eddy viscosity as a function of stability and vertical shear of the 
v/ind« 
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The Coefficient of Eddy Viscosity as a Function of Height and Time 0 
The horizontal equations of motion for turbulent flow where the eddy 



stresses are assumed to be a function of height and time only and the 
molecular or viscous stress is considered negligible are 

h- + V'V« = f(v-v s ) + JL 

c>t r ^ 

Here "T^x; and are the components of eddy stressj 

and 'V~q the components of the geostrophic wind} , the mean 

wind velocity vector of which U andiT are the horizontal components; 

Y y the density; and 4~ , the coriolis parameter « 

The eddy stresses are assumed to be related to the mean wind shear 
by the expressions 






J t)Z. 



( 2 ) 



The geostrophic departure method employs observations of the mean 
velocity and geostrophic wind to evaluate all the terms of eq (1), so 

tTzx 



as to yield the distribution of 






and hw 



This method 



gives the distribution of stress and eddy viscosity provided the wind 

bu , 

components have a maximum at some levels h x and hy so that at h x ,xr~-U 
and at h. r ^ — - - - 0 . 

0 Z- 

By multiplying eq (1) by ^ and integrating from h x to an 
arbitrary height z, the distribution ofi^* is obtained. Similarly 
the distribution o f 1 Z<j is found o Then eq.(2) leads to two independent 
estimates of the distribution of eddy viscosity 0 



( 5 ) 





Assuming no discrepancies in the procedure or data,, 0. 




Blackadar ca used the geostrophic departure method on a 



composite wind distribution obtained by averaging the winds at two 
stations (Oklahoma City, Oklahoma and Wichita, Kansas) over 2k 
selected, not necessarily consecutive, one-day periods 0 Figures 1 



periodo Although great care was taken in determining each of the 



’’the independent estimates of and did not 
agree and negative values of each were numerous and 
frequently large D It is assumed that the reason for 
the poor performance is in the inaccuracy of the 
geostrophic wind 0 n 

For this reason, he devised an indirect method for obtaining the 
values of and-^Vg which would minimize the sum of the absolute 
magnitudes of at the four levels ( 950 * 900* 850 and 800 mb) 0 

He forced JA K to be equal to at a single level and at the two 
standard radiosonde times J then the geostrophic wind components were 
computed*, Next utilizing the thermal wind relationship and assuming 
the temperature sounding to be accurate, the values of Ug and^ were 
found in the vertical e The variation of geostrophic wind at times 
other than the standard radiosonde times was determined on the 



and 2, reproduced from Blackadar 5 s D1 data, show the eastward 
and northward components of the average winds during the 2k day 



terms of eq (l), 




pointed out that 



k 




(Sd3J.3W) -LH9I3H 
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The eddy 



-raises and eddy 



3 n r e o onp ut « ” u t - ‘-0 .•/«* * e r he i ght 



.ntervals and three how .Vumals, The values of and ^ v h- 

so obtained were generally x -gr r. w.. ent at. all t mes 0 Blackadar 1^1 
aeci ded 



Fig 0 5 shows Blackodar ] s .1 computed values of as a function of 
height and tine. 

It is desired to find a function (z,,t) which will yield these 
values of /i 0 Due to the nerieenc nature of the solutions, fig* 5 
must be modified in order to show the sane values of t ? for 
corresponding heights at the beginning and end of the 24 hour periodo 
The alteration of -fig* 5 was minimized insofar as possible, and the 
modified figo 5 is shown as fig 0 r K 

The following approach to the problem was used. Determine the 
time series which fits the observed values of , for 12 times 
( 2 hour intervals) at each of the given heights 0 y 200 >. ^00, 6OO0000 
i£00 meters* An example for the height of 200 meters is 



Here sine and cosine terms belonging to the same harmonic have been 
combined^ Aq( 2) is the mean value* Aj (2) is the amplitude of the 
first harmonic* T^(2) is the phase angle of the first harmonic in 
degrees and t is time in hours after t ~ 0 o (1500 CST}* expressed in 



n 



to dispose of the lemaimag differences* the best 
imate - was assured to be given by 



estimate 




2 - (2,t) = A c (2) * A- ; (2) cos [ (t ♦ Tj (2) 1 

r. . , , ^ . . - (^) 
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a \v\tW 



degrees How, having xO equations of the same form as eq (4), fit a 



continuous function of z to the 10 values of A and T for each harmonic 



The harmonic analysis at each even 200 meter height was accomplished 
yielding the values in Table 1. 

Fig, 5 shov/s values of Aq (the mean), A] and Ag versus height; 
fig, 6 shows Aq, and Ac- versus height, Fige 7 shows 7^ and Tg 
versus height and fig, 8 shows Tj y and Tq versus height, 

A diagonal difference table was computed for Aq values versus 

g 

height, The second differences of A Aq^ showed the most nearly 
constant values. Although the second differences were not very close 
(ranging from to 128), a parabola was fitted to the values of Aq 
by the least squares method. As expected, the parabola did not fit the 
data to within sufficient accuracy. 

Next it was noted from figs, 5 and 6 that a possible orderly 
decrease in the ratios of A^/Aq existed at each height ax i increased 
from 1 to 5° A plot of the ratios A^/Aq versus height was drawn 
(see fig, 9) o Utilizing these ratios the equation foryM(z,t) would 
take on the form 



which makes a significant contribution. This would yield the final 



form of the equation 




( 5 ) 
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TABLE 1. 



Results of harmonic analysis at 10 levels 0 



Z is height above surface (m) , the mean coefficient of eddy viscosity 
(g cm~^sec~^) , the amplitude of the i ^ harmonic (g cm" 1 sec ^ ) and 
T x is phase angle of i ^ harmonic degrees 0 ' ' 



z 


A 0 


A 1 


Ti 


A 2 


t 2 


a 5 


t 5 


a 4 


t 4 


a 5 


t 5 


a 6 


0 


204.1 


255.4 


26.1 


86.2 


40.2 


4o.o 


58.7 


22.9 


54.9 


7.6 


96.8 


2.6 


200 


118.5 


158.8 


19.7 


74.8 


54.9 


55.6 


61.2 


24.4 


85.7 


19.4 


125.8 


8.8 


4oo 


75.8 


87.8 


12.9 


54.5 


51.9 


24.8 


69.1 


16.6 


115.5 


18. 5 


156.0 


9.5 


6oo 


87.0 


105.5 


15.7 


50.0 


54.5 


19.4 


69.4 


15.5 


125.8 


25.7 


165.1 


12.2 


800 


144.5 


191.9 


11.9 


98.7 


27.6 


45.8 


45.7 


18. 5 


106.5 


26.0 


152.2 


14.4 


1000 


266 


577.6 


7.7 


255.5 


16.0 


158.1 


18.6 


95.0 


19.5 


68.6 


18.7 


51.2 


1200 


461.2 


690.0 


5.5 


471.5 


8.2 


551.1 


8.5 


248.4 


6.5 


196.6 


6.9 


95.7 


1400 


758 


1162 


2.2 


806 


7.1 


560 


8.5 


560 


8.0 


258 


8.0 


114.6 


1600 


lloO 


1761 


2.1 


1169 


8.6 


756 


11.8 


458 


15.6 


266 


15.2 


105 


1800 


1690 


2518 


5.4 


1620 


18.1 


984 


51.0 


559 


44.0 


249 


45.5 


60 



( 1 ) 



All surface values are the result of extrapolating Blackadar*s ■ 1 
curves down to the surface 0 
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Her 3 Aq(z) is tne function of height needed to fit the observed values 
of the mean yM at each height ~ ®i(z) is the function of height needed 

to fit the ratio A x /Ao at each height and T x (z) is the function of 
height needed to fit the observed values at each height 0 If the 
curves of fig. 9 could be fitted by a relatively simple function of 
height vice an estimated ^ih or 6th degree polynomial^ there existed a 
possibility of a series analysis and possible further simplification 0 
One noteworthy aspect of this approach should be observedo It appears 
that a fairly simple function of height (a e g 0} , a parabola) could be 
fitted to the ratio curves (fig 0 9) with reasonable accuracy if a 
layerized approach is followed 0 The dashed lines on fig* 9 represent 
the probable curves which would fit a 2 layer model e The division of 
layers appears to be in the vicinity of 600 to 700 meters 0 This 
conceivably could be of physical significance * 

Next it was noted that Aq (fig© 5) increases non-linearly with 
height above 460 meters 0 Further, the point z = 450 meters , for which 
A © = JJg.cm sec 9 appears to be a point of syrnmetry 0 The axis was 
translated to this point and Aq versus height was plotted on log- log 
paper; this plot coincided almost exactly with a straight line c 
Ey the least squares method^ Aq(z) ~ az n * b was determined 0 
The equation for Aq was 

A Q (z) - 6„c |z - 4.6 | 2a 4 75 (7) 

where z is height above the ground in hundreds of meters 0 Table 2 
shows a comparison of the computed values of Aq and observed values 
(from the harmonic analysis) 0 
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TABLE 2. 

Aq( 2 )(S cm ^sec - ^) computed versus Aq(z) observed. 



HEIGHT 

Z * x 10 -2 

meters 


COMPUTED 

A 0 (Z) 


OBSERVED 

A 0 (Z) 


PERCENT 

ERROR 


0 


225 


204. 


10.2 


2 


118. 


118. 


0.5 


4 


75- 


75- 


5.6 


6 


85- 


87 


2.0 


8 


155- 


-3" 

i — i 


6.5 


10 


286. 


266 


7.6 


12 


487- 


461. 


5.7 


14 


759- 


758 


0.2 


16 


1105 


ll6o 


4.7 


18 


1524 


1690 


9.8 
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Sinilarly A^(z) a ^d A 0 (z) were obtained. A comparison of A^(z) 
computed and observed is shown in Table and for A 0 (z) in Table 4. 
The equations developed by the least squares method are as follows: 



A j (z) =8.4 
A 2 (z) = 2.5 



z - 4.5 | 2,1 5 + 81 


(6) 


z - 4.0 | + 55 


(9) 



When locating the approximate point of symmetry of the Ag versus 
height curve (fig. 5) » it was necessary to employ the method of least 
squares to determine which of three logically chosen positions for the 
height axis yielded the least difference between A 0 (z) calculated and 
A^(z) observed. In spite of this effort, relatively large differences 
were observed in the range, surface to 400 meters (see table 4), where 
the symmetry was weakest; however, the fit appears to e good above 
oOC meters. 

A linear approximation was used to fit T^(z) and T 0 (z) (see fig. 7) 
The equations determined by the least squares method are 



T-j (z) =-1.20 z + 21.4 
T 0 (z) = - 1.81 z + 59 

C- 



( 10 ) 



where z = height above surface in hundreds of jaeters and T-^(z) and ? 0 (z) 
are in degrees. Table 5 shows a comparison of the phase angles computed 
by eq (10) for the first and second harmonics and the observed phase 
angles . 
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TABLE 5- 



A : (z)(g :m*"^sec”^) computed versus A^(z) observed 



HEIGHT 

^meters x ^ 


COMPUTED 

4p(Z) 


OBSERVED 

A : (Z) 


'DTp'D OT? rT ’'p 
rJJxWOA'ii 

ERROR 


0 


295- 


255. 


25.9 


2 


i — 1 

- 4 * 

i — 1 


15 1 


1.7 


4 


81. 


8$ 


7.5 


6 


165 


105 . 


0.5 


8 


205 


194 


6.8 


10 


408. 


577. 


8.2 


12 


719. 


690 


4.5 


14 


1142 


1162 


1.7 


16 


1685 


1761 


4.5 


18 


2540 


2518 


7.1 
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TABLE 4. 



A 0 (z)(g om - ^sec - ^) computed versus A 2 (z) observed 



HEIGHT 

^meters x ^ 


COMPUTED 
^2 (Z) 


OBSERVED 

A X (Z) 


PERCENT 

ERROR 


0 


110 . 


86 . 


21.6 


2 


48. 


74. 


5^.9 


4 


55- 


54. 


55.5 


6 


48. 


50. 


2.6 


8 


110 . 


98. 


11.8 


10 


250 


255. 


2.4 


12 


4 SO. 


471. 


4.8 


14 




806 


6.8 


16 


1157 


1169 


1.0 


18 


1674 


1620 


, 5 -> 
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TABLE 5. 



T^(z) coaputed versus T^(z) observed 
and Tg(z) coaputed versus TgCz) observed for 10 levels. 
Tj(z) and TgCz) in degrees. 



^aeters " 


T]_ (z)obs. 


T]_ (z)coap. 


T2' , .z)obs. 


(z) comp. 


0 


26 . 


21. 


4 o . 


59 


2 


19 . 


18 . 


IS 


55. 


4 


13 


16 . 5 ; 


1%, 


IX 


6 


IH- 


14. 1 


7 A. 


28. 


8 


i% 


12 . 


27 . 


24. 


10 


7 . 


9. 


16 . 


20. 


12 


5- 


7 


8. 


17 . 


14 


2. 


4 . 


7 . 


15 - 


16 


2. 


2.0 


8. 


! 

10 . 

i 


18 


5. 


- 0.4 


18 . 


6. 
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Upon examination of the remaining amplitude coefficients and their 
corresponding phase angles (see figs. 6 and 8), it appeared that the 
functions of height needed to fit these values to the desired degree of 
accuracy would he relatively long and too complex to warrant further 
effort in this phase. 

Also one of the objectives v. T as to see hov; complicated a function of 
height and time it would take to yield the observed values of in fig* 5* 
The equation foryM (z,t) thus far obtained is 




where z = height above ground in meters x 10 “ and t is time in hours 
after t = 0 (1^00 G3T) . Values ofy^f (z,t) were computed using the 
Ilarchant desk calculator for time intervals of one hour and every 100 
meters height. Fig. 10 shows computed values of /"/plotted on a 
height-time coordinate system. 

Comparison of fig* 10 and the desirod^M distribution from fig* 4, 
shows fig* 10 to contain a fairly representative set of isolines for 
from 1000 C3T to 2000 C3T from the surface to 1800 meters. Both 
figures show the largest order of magnitude for the coefficient of 
eddy viscosity at approximately 1^00 C3T at all levels. In addition, 
both figures show a minimum f or /A in the vertical at approximately 
400 meters between 1000 GST and 2000 C3T withyM decreasing from the 
surface to 400 meters throughout the period and with ft increasing 
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rapidly with height above 400 meters from 1200 CST to 1300 C3T. In 
fig. 10, however, the decrease of /A with time at any level above 1000 
meters is too great in the vicinity of 1900 C3T, as is the increase of 
fA with time too large at any level above 1000 meters in the vicinity 
of 1000 CST. Also, the vertical gradient of /A above 4C0 meters at 
l^OC GOT is less than that observed on fig. 4. The latter could be 
rectified by computing the remaining harmonics which would tend to 
increase the values of fi computed at 1^00 CST (t = 0) so that they 
more closely approach corresponding values of observedyM (fig. 10). 

Between 2000 CST and 1000 CST, surface to 1800 meters, a 
comparison of fig. 4 and fig. 10 shows fig. 10 to have a completely 
unacceptable representation of the distribution ofyAf . Fig. 10 shows 
a column of minimum^ values extending vertically from the surface to 
1800 meters at about 2200 CST, with no phase lag, and another minimum 
column, extending from the surface, at 0^00 CST, to 1800 meters, where 
the minimum occurs at about 0700 CST. These minima are separated by 
a relatively strong maximum, surface to 1800 meters, occurring at 
0200 C3T. In contrast, the desired form, fig. 4 , shows a strong 
single minimum extending from the surface, where it is estimated to 
occur at 2pOC CST, to 1800 meters where it is observed at approximately 
0500 CST. 

The double maxima, double minima vice the desired single maximum, 
single minimum, are attributed primarily to the fact that the mean and 
only two of the six harmonics were utilized in deriving the values of 
fA for fig. 10. Fig. 11 is introduced to lend credence to the above 
remark. This figure shows the resultant of the first two harmonics 
versus time for height equal to 1200 meters. The resultant 
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distribution utilizing the mean and the first two narnomcs as 
determined by the harmonic analysis, not from fitted values, tends to 
show the following characteristicsSan absolute maximum at 1500 C3T, a 
relative maximum at 0250 C3T, and two minima, one at 0700 03T and the 
other at 2200 C3T 0 

Each of the above maxima and minima are observed on fig. 10 at 
the corresponding times and the magnitude comparison is as follows; 

TIME jUi from fig. 10 Ja from fig, 11 

1500 C3T I650 1620 

2200 C3T -85 -l4o 

0250 C3T 215 240 

0700 CST -124 -140 

Similar analysis at other heights would yield similar results 0 

This would tend to show then that the fitted values have not 
distorted the picture of distribution with time and height except to 
a slight degree and that the cause of the deviation of the computed JA 
values (fig. 10) f ron Blackadar * s no observed values (altered 
fig. 4) must be the remaining four harmonics 0 Further,, when using 
eq (11) and including the remaining four harmonics, the values of 
computed compare very favorably with the observed values 0 
Pertinent to the variation of with height, Sstoque ^5*1 
has recently approximated the distribution of JA by a linear decrease 
of above a maximum value found at a height of about 5 0 meters. 

He proposes a two layer model for the frictional layer, the lower 
layer, called the surface sub-layer, extends to approximately 50 
meters, and an upper layer, called the transition zone, extends to 
approximately 2000 meters. Utilizing this JA distribution, he has 
numerically integrated the equations of notion and heet conduction 
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equation obtain the 24 hour distribution of wind and temperature • 
The diurnal wind distribution thus obtained by Zstoaue appears 
to be in fairly good agreement with the diurnal wind characteristics 



to 



data 



observed on several occasio ns c However, for Blackadar's 
(see fig 0 p)$ it vas necessary to attain an average height of 
approximately 11^0 meters before Ja varied linearly (but increased) 
v/ith height to within a 10 percent error of the observed value of JA 0 
Further, the slopes of the lines representing the linear variation of 
JA v/ith height (for different times) varied up to 55° • 



28 



as a Function of Stability and 



The Coefficient of Sddy /iscosity 
Vertical Shear of the Wind* 

Wind data for this phase were obtained from figs* 1 and 2, 
reproduced from Blackadar lf ll * Temperature data for the 24 selected 
days were obtained from photo copies of WBAII-55 for the months of June, 
July, August, and September 1952 for Wichita, Kansas and Oklahoma City, 
Oklahoma, which were furnished by the National Weather Records Center, 
Ashville, North Carolina* 

All temperatures at the reported levels (surface, 950, 900, 850, 
800 and 750 mb) at both cities and at both sounding times (1500 GCT 
and 0p00 GCT) were converted to potential temperature* The surface 
potential temperature was obtained by the Arowagram and upper level 
potential temperature was obtained from the Smithsonian Meteorological 
Tables* 

For the 0p00 CCT sounding and at each level, the average potential 
temperature was obtained by adding the 24 observations at V/ichita to 
the 24 observations at Oklahoma City and dividing by 48* This yielded 
the average sounding for the two cities over the 24 day period at OpOO 
GCT* A similar procedure was followed to obtain the average height of 
the 950, 900, 850, 800 and 750 ~b levels. All heights are in terns of 
0*98 geodynamic meters* The same procedure was followed at 1500 GCT* 

The average heights and potential temperatures obtained for the 
two times are shown in table 6* Next the stability factor j 

J_ bo „ \ 0 * 

e bz 

-,.' 2.3 computed for the center of the averaged layers using finite 
differences for the partial derivative of & . Stability factor 
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TABL2 6 



Twenty four day average heights (in 0*98 geodynamic meters) 
and potential temperatures (degrees K) for the reported 
levels of concern. 



Level 


Z 

meters 


e 

degrees K 




Level 


J. yUW UW 

Z 

meters 


e 

degrees 


Surface 


0 


5 00. g 




Surface 


0 


501.8 


950mb 


188 


505.0 




950nb 


201 


502.5 


900 mb 


^59.1 


506.3 




900nb 


672.5 


505.1 


SyOnb 


1156.5 


507. j» 




850nb 


1168.2 


507. C 


80Cmb 


1677 


509. 3 




800mb 


1687.9 


508. (0 


790mb 


2228.2 


511.5 




750nb 


2257.6 


510.7 



VN 



versus height of center of layer for the two tines is shown in table 7o 

>v 

Values of the components of the vertical vand shear, a.nd _ , 



* bz- “““ b£ 



for each 100 meter level were then obtained from fijto. 1 and 2 
utilizing finite differences for the partial derivatives „ The shear 



faotor '0M + ©K 

level for the two tines , 



was computed next for each 100 meter 



Graphs of^ < observed versus the stability factor and of/M observed 
versus the shear factor were then constructed in an attempt to find a 

correlation between these factors andy^ « One apparent correlation was 

% 

that in the lower layers (below 600 meters) where pronounced stability 
variations occurred, strong stability was associated with a small , 

and vice versa* Above about 600 meters the contrast in stability factor 
was considered negligible* Moreover, there was insufficient difference 
between the ft values at each level for these two times (OfOO C3T and 
2100 C3T) , nor between the shear factors in the upper levels (above 
600 meters), to detect any correlation between wind shear and the 
values ft observed . In the lower layers, when the shear appeared to 
be large, the value of yM was also large* In order to see if a 
correlation between wind shear andy^ existed in the upper levels, it 
was necessary to examine the data at tines- of minimum and maximumyM 
(ioCo 0p00 C3T and 1^00 GST)* Computation of the required shear factors 
showed a strong correlation between low shear values and highyW values 
and vice versa for all levels* This did not contradict previous results 
above 600 meters but it was in opposition to the relationship found 
between shear and/< below o00 meters at 0900 and 2100 GST* 

Since an expression throughout the vertical for ft in terms of 
stability and wind shear is a most desirable form, a compromise appeared 
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TABLE 7. 



Stability factor , versus height of center of layer for 

© c> £ 

24 day average soundings at OyOO GOT and 1500 GOT. 





0500 GOT 






1500 GCT 




Layer 


Height of 
Gen 4, 

Vmo-oers ) 


J-dS.lO 4 

(meters ~ { ) 


Layer 


Height of 
center of 
layer 
(meters) 


1*6.1 
( meter. 


Surface to 950ml 


94 


. 7498 


Surface to 9^>0mb 


100 


.0807 


950nb to 900mt 


424 


.0826 


950nb to 900rab 


457 


.1949 


900mb to 850ml 


908 


.0865 


900mb to 850mb 


920 


.1265 


850mb to 800mb 


1417 • 


.1065 


850nb to 800nb 


1428 


.0969 


800mb to 750mb 


1955 


rH 

KA 

rH 


800nb to 7 50® h 


1965 


.1250 



.10 
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to be necessary, i*e 0 , a two layer model with yH a function of the 
stability factor below 600 meters and/^f a function of the shear factor 
above 600 meters. Now, from visual inspection of figs. 1 and 2 at 
1500 GST and from prior shear computations the shear can be seen to 
be very small and practically constant throughout the layer surface 
to 1800 meters; yet the value ofytf (see fig. 5) increases rapidly 
from 600 meters upwards. Therefore, a second compromise from the 
originally intended solution appeared to be in order, i.e., above 
600 meters let yW vary with the shear factor and with height. 

Further, due to the lack of temperature data below 600 meters, it was 
decided to utilize the mean (z) with respect to time in an attempt 
to more accurately describe the yH observed values and also perhaps to 
extend the coverage below oOO meters to the entire 24 hour period 
rather than merely the times of the two temperature soundings. 

The general approach used in the solution of the problem is as 
follows: 

1) Assume (for both layers) 



where is the observed value of /M for any point in time and space, 



is the difference between the observed value ofy^and the mean value 
of yM at that height. 

2) Find the best function of height which yields the /M (z) 
distribution. 




( 12 ) 



yH (z) is the mean value of all observed at the level z, and &/A 



5) Account for A/^ , the deviation from the mean, 
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a) below oOO meters by a function of stability factor, 

b) above oOO meters by a function of shear factor * 



A function was found in Section 2 to represent the mean^Af at any 
height (see eq (7)); however^ this was based upon values obtained from '"V 
the isolines ofyH in fig* 4. For this section an attempt was being 
made to find the function which yielded the computed values plotted 
on fig, 5> inasmuch as they were in part computed from the vertical 
shear of the horizontal wind. A similar approach to that used in 
deriving eq (7) was in order and yielded the following foryM (z) 



where z is height above the ground in hundreds of meters. Due to the 
relatively good fit, this equation was taken to apply in both layers. 



Four of the points plotted are from observed data for the period 
(table 7) and the fifth point is an estimated value for a height of 
150 meters at 1500 C3T. This value was derived from 12 soundings 



^ (z) ft. 5.6 |z -4.8 | 2,2 + 65 



(15) 




computed with eq ( 15 ) from the values of ^ observed on fig. 5* 

versus stability factor is shown as fig. 12. 




taken during the Great Plains Turbulence Field Program 




August and September 1955 and under similar atmospheric flow 



conditions. The stability factor obtained was .0055 m \ and the 
A/f corresponding to 1500 C3T and 150 meters is 7 4f>g cm~^sec~^ 0 
The curve appearing on the figure has been faired in and does not 



represent the fitted function given below. 



Of various functions of the stability factor tried, that 
which yielded the smallest "least squares" error for was the 

hyperbola, 
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35 



- 1 21 



5 



(3* ) + i [ 200 , 5 ^ - 7^00 (3 0 ) + .12J (3 e ) 2 J 



where (3. ) --L }~S-. x 



A • 

^ - 



(14) 



Combining eq (15) and eq (14) then yields the desired form for 
computing/^ , i„e, 



yH computed - 5-60 



z ~ 4„8 



2 ° 20 + 186 - 5«6 ( 3 & ) 



— j 1 

200 ,500 - 7400(S & ) + 125 {3 & ) 2 J 



( 15 ) 



Here z is in meters x 10 

Table 8 shows a comparison of yM computed and /^observed for tho 
layer below 600 meters 0 The largest error occurs at 2100 C5T and a 
height of 424 meters 0 It should be noted (see figs, 1 , 2 ) that at this 
point the shear is relatively strong which, according to eq (3 X would 
tend to yield a small value for /Y 0 A similar analysis applies for the 
next largest error, percent wise, at 0900 C3T and 457 meters 0 Here, 
where the. observed value is greater than the computed value, a low 
value of wind shear is observed 0 From the preceding analysis, it 
appears that shear is becoming a factor in the upper levels of the 
lower layer o In this transition zone from approximately 400 to 600 
meters, a better representation of might be obtained by considering 
shear as well as stability*. 

For the layer above 600 meters, a graph was drawn showing height 
versus shear factor and the associated value (next to each point) 
needed to yield the observed value of for the four times shear 
computations were made , This is shown as fig 0 1% 
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TABLE 8 



Comparison of fA computed and y •' observed(g onf^sec - ^), 

surface to 600 meters 0 



Tine Height above Surface 

’C3T) (meters) 


M observed 
' - 1 _ 1 

g on 3ec x 


/ ' computed 
g ca*"^sec"*^ 


C\ 

o 
o 
< — 1 
OJ 


60 


60 


2100 424 


27 


4 9 


0900 100 


l4o 


155 


0900 4p7 


* 

24 


16 


1500 150 


480 


45o 
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As previously stator. j.u -./as desired to account for AyM values as 
a function of height and shear* The approach to the problem was 
similar to that used in section 2* An equilateral hyperbola was fitted 
to the observed values of shear and ts.fi at each level (oOO, 800, 1000 
o.odSOO meters) to provide the A/^ distribution as a function of 
shear factor* The result is 



A/*. 



Cz 



+ b: 



Z - 



IfeNlSr 



(16) 



where t\j^ z is the difference between the observed yM and yM (z) at 
height z above the ground; c z and b z are the constants determined by 
the least squares method at the z level* Fig* 14 shows c z and b z 
versus height* The curves appearing on the figure have been faired 
in and do not represent the fitted function given below* 

IText the c z and b z values deteimined at each height were fitted 
with a continuous function of height* The functions of height 
representing each of these is 



c(z) = *22 (z) fc 

b(z) = *10 z)^ 0 ^ 

*2 

where z is height above ground in meters x 10 • 

A comparison of the c z and b z values observed and the c(z) and 
values computed is shown in table <?<= The resulting function of 
and height representing the values at 600 meters and above 



(17) 

(18) 



b(z) 

shear 



is 





( s ) 
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( 19 ) 




4o 



The constants 



c„ and b z versus c(z) and b(z) at each level 



z 

~2 

meters xlO 


c z 

g crn"^sec“2 


c(z) 

-1 _3 

g cn sec ' 


b z 

-1 -1 
f zrz Gee 


b(z) 

g cnf 1 sec~ 1 


s 


514 


509 


-56 


-51 


3 


625- 


982 


-119 


-]4o 


10 


5604 


2440 


-520 


-502 


12 


8596 


4980 


-609 


-574 


14 


7457 


9531 


-870 


-970 


16 


14,572 


16,115 


-1590 


-1541 


18 


24,147 


25,840 


-2578 


-2518 



4l 



<h) u? 



and the function representing M in the upper layer is 



yM = 5 <=6 (z - 4 C 8) 2 ” 2 



65 + c2; 




( 20 ) 



A comparison of the computed Rvalues and the observed values is 
shown in Table 10* 

The distribution of fA computed indicated in Table 10 compares 
fairly well with the observed values* The salient features of the 
observed^ distribution are preserved by the computed pattern, i 0 e 0 , 
a maximum at each level at time 1^00 C3T (one exception at 600 meters) 
and a minimum at each level at time 0J00 CST. In most instances the 
order of magnitude of the observed and computed values are the same 
(five exceptions out of 28 observations). The value observed at 
0900 GST, height 600 meters, is approximately one fifth that computed 
by shear and height function* This discrepancy may partially be 
accounted for by the higher stability factor observed in the vicinity 
(see table 7)® The negative values computed for heights 1200 and 1400 
meters at 0J00 C3T are primarily a result of the curve fitting process 0 
(3ee fig. 14 and table 9)* 
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TABLE 10 



Comparison of computed coefficient of eddy viscosity ,/f c (g cn” 1 sec~^) 
and observed, u Q for four times in upper layer , 



Height 

^meters) 


Me 


(>o) 


Me 


(A 0) 


A 




A 


fa) 


1000 


7607 


(7020) 


455 


(480) 


105 


( 157 ) 


1100 


(528) 


loCO 


6709 


(5900) 


734 


(552) 


56 


( 97 ) 


509 


(208) 


l4oo 


5800 


(4680) 


142 


(200) 


-29 


( 85 ) 


177 


(155) 


1200 


1811 


( 3000 ) 


44 


(H5) 


-7 


(57) 


120 


(112) 


1000 


836 


( 1230 ) 


4l 


(83) 


7 


(20) 


75 


(81) 


800 


577 


(588) 


24 


(^5) 


15 


(17) 


46 


( 80 ) 


600 


99 


(142) 


48 


(58) 


26 


(28) 


228 


(42) 



21Z OJZ 0$Z l, r Z 

pro 21C OJC OfO 






4, Conclusions * 



From the approach followed in this investigation a simple 
representation of the diurnal variation of the coefficient of eddy- 
viscosity with height was not apparent * 

It has been shov.m that the function of height and time 
representing the mean and first two harmonics is inadequate to show the 
salient features of the observed distribution. 

In order to represent the observed distribution of the coefficient 
of eddy viscosity by a function of height and time, it would take a 
relatively long series of terms, estimated to encompass at least five 
harmonics, involving exponential functions of height, fourth or fifth 
degree polynomial functions of height, and trigonometric functions of 
angles dependent on time and functions of height, 

A two layer model is developed in an attempt to account for the 
observed distribution of the coefficient of eddy viscosity as a 
function of stability and shear of the wind. It is found useful in 
the lower layer, and essential in the upper layer, to include height 
above the ground as one of the parameters determining the distribution,, 
The variance of the coefficient with stability in the lower layer 
appears well founded, as does that bet ween shear and the eddy 
viscosity coefficient in the upper layer; consequently, a fair 
representation of the general characteristics of the observed diurnal 
variation of the coefficient of eddy viscosity is obtained from the 
empirical formulae. representing the two layers 0 

Further investigation utilizing more frequent and detailed wind 
and temperature observations would bo likely to yield better results 
and might eliminate the need for the two layer approach,, 
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